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ABSTRACT 


Giaever and Zeller (1968) have shown that small 
Metallic partrelessin che barrver of a thin’ f21m 
tunnel junction cause a strong conductance dip to 
occur at zero Dias in the tunneling characteristics. 
The characteristics are modified slightly if the par- 
ticles are superconducting. These effects have been 
studied for Sn and Pb particles in the barrier of 
A1-I-Al tunnel junctions. The values BE theytransi— 
tions temperatures of Sn particles are found to be 
adhe, ha ROA OB yes. Ziel a, 0D) ~ 07K LOL Particles of average 
thicknesses 65 & and 130 A respectively. The difference 
between the values of the transition temperatures of Pb 
particles of average thicknesses 60 A, 80 & and 130 A 
CouLG snot be" distinguished. ~lt was found to be 7. 44.05°K 
for all the three sets. The values of the energy gaps 
aLeeLOund «co .Dew2.59eMVge2 05 MV, sanGe2./1 my vespectively 
aol O082 Kes toewpartilcie Size diStrebucion for vay set iof 
Pb particles Pes en from electron micrographs in 
order to test the Giaever-Zeller theory. Small effects 
due to the phonon spectrum of Pb particles were also 


observed. 
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CHAPTER 1 


INTRODUCTION 


Since Giaever's (1960) work on the experimental 
aspect of tunneling, there has been a sequence of 
measurements on both superconducting and normal metal- 
insulator-metal junctions. The results have continually 
puzzled both theoreticians and experimentalists leading 
them into one problem after another. 

Minimum oe for the occurrence of superconductivity 
has interested many workers. Anderson (1959) in his 
"Theory of Dirty Superconductors" predicted a lower size 
limit for the occurrence of superconductivity when the 
Spacing of the electron eee levels is less than the 
energy gap. In his discussion of "Electrostatic Effect 
and souperconductivity in Tiny Particles» Markowitz (1967) 
modified the BCS Hamiltonian by adding on a Coulomb term 
and then minimised the ground state energy by considering, 
among other parameters, the sample size. These considera- 
tions led him to estimate the minimum size of a supercon- 
ductor. Parmenter (1968) has derived an expression 
relating the dependence of the energy gap and transition 
temperature upon particle size. His calculations lead to 


Gather too high values of critical size. 
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Zeller and Giaever (1968) have reported experimen- 
tal results on small superconductors and proposed a 
capacitor model to account for the observed zero-bias 
resistance anomalies. They also observed particles as 
small as about 25 A superconducting. 

This thesis contains a report of the energy gaps 
and transition temperatures of small superconducting 
particles and some differences between direct tunneling 
and tunneling via particles in both normal and super- 
conducting states. The size effect in superconductivity 
manifests itself in the energy gap and transition tem- 


perature which are different from those of a bulk material. 
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CHAPTER 2 


MICROSCOPIC THEORIES OF SUPERCONDUCTIVITY 


The BCS mnSOLy. 


In the course of measurements of the decrease of 
electrical resistance of metals when they are cooled, 
Kameriangn Onnes inelotl Lound that mercury lost al. or 
its resistance suddenly at about 4°K. Onnes called this 
state of no resistance the "superconducting state". For 
many years this phenomenon defied theoretical explana- 
tions. The first attempt was made by Gorter and Casimir 
(1934). They suggested a two-fluid model which stipulates 
that a fraction x of the metallic electrons remain normal 
at the transition temperature, Tor while the remainder 
(l- x) is condensed into the superconducting state. 


Recording sto “this model, the traction x 1s given by. 


the specific heat in the superconducting state is 


dh 
(© = Shp oe Die, 
m2 

Cc 


where y is the Sommerfeld constant, and the critical 


field is given by 
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These results were not in perfect agreement with experi- 
mental results: 

The next step was made by London (1935) who derived 
"London equations" by introducing a characteristic cons- 
tant Ay, called the penetration depth. But the Meissner 
effect required that a field already inside the material 
in the normal state should be expelled when the material 
becomes superconducting and this was not accounted for. 


London added another equation by postulating that 


2.4 


London's results together with those of Gorter and Casimir 

were used to fit experimental results but not successfully. 
With the discovery of the isotope effect by Maxwell 

(1950) and Reynolds et al (1950) showing that lattice 

vibrations were involved in the superconducting state, 

a new approach was needed. Pippard (1953) introduced 

a non-local theory as a modification of the London theory, 

with a new parameter” € called the coherence length. In 


terms of the microscopic theory, § may be considered as 
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the typical size of the "Cooper pairs". In the London 
equation 
d(x) =- —S5 A(t) =- = air) , 2.6 
a = cm 
ATA 
L 
where 
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Pippard introduced a dependence of the penetration depth 
on the electronic mean free path, L, such that 
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where E5 is a-constant of the superconductor in question. 
Derivation for the field penetration into a semi-infinite 


slab gives the field 


H(x) =" He e X/4 


where 


Earlier on Pippard had shown that under low temperatures 
and high frequencies the mean free path may be longer 
than the skin depth, the electrons become less effective 


as current carriers and Ohm's law 
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is no longer a reasonable approximation, but rather the 
current at a point must be determined by the total field 
over the distances of the order of the mean free path. 
Reuter and Sondheimer (1948) derived an expression for 


this anomalous skin effect of the form 
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vg 

As analogous to equation (2.10), Pippard proposed a rela- 
tion for the electromagnetic response of a superconductor 


of the form 
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A more successful explanation for the superconduc- 
ting phenomenon has been presented by Bardeen, Cooper and 
Schrieffer (1957) known as the BCS theory. The BCS theory 
was set up to show that superconductivity was bound up 
with the behaviour of the conduction electrons; while the 
DbEoperties OL ther crystal e@lattice doesnot change, athe 
GoyV Stale lattice plays an important; party and, that the 
superconducting-to-normal transition is a second order 
phase change. 

Frohlich (1950) had postulated that an electron- 
phonon interaction was able to couple two electrons such 


that they could behave as if it was an electron-electron 
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interaction. In this picture, an electron with momentum 
Ki interacts with the lattice emitting a phonon of momen- 


tum q and leaving the electron scattered with momentum ki- 
The phonon is then absorbed by a second electron with 


2 


of the phonon, figure 2.1. By the law of conservation 


momentum k, which changes to momentum k,. on absorption 


of momentum, 


kK, =k, +@ 
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While the momentum and the energy are conserved between 
the initial and the final states, anything can happen in 
between the emission and the absorption of the phonon. 
If the life time, At, of this intermediate state is very 
short, there will be a large uncertainty in its energy 


AE according to the uncertainty principle 


AE.At = h 


in which case the energy need not be conserved in the 
emission and absorption process. Cooper (1956) carried 
forward this idea saying that there are correlations 
between electrons of opposite spin and these correlations 


are produced by the interaction between the electrons. 
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While in a classical model a very strong two-particle 
repulsion between two electrons implies that the two 
electrons would ineverggetyvery close to one another, at 
very low temperatures the system can gain energy by 
going into a highly correlated state, Cooper (1959). 
Cooperspost! aktedgthek, thesinteractromipreducing the 
correlations is between electrons via the lattice. 

When an electron collides with a lattice wave, 
it may be scattered producing some resistance. At 
T = O0O°K, there are no lattice vibrations and ithe phonon 
part of the resistance goes to zero. But an electron 
can excite a virtual lattice wave which in turn inter- 
acts with another electron. Alternatively one can say 
that when an electron moves through a lattice of in- 
finitely heavy lattice points, the lattice is polarized 
and another electron moving through it at a later time 
feels a distorted potential due to the polarization. 
The distorted potential is a function of position and 
time and shows up as a retarded interaction between two 
electrons. 
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correlation produced by the interaction between electrons, 
r is the relative coordinate of the two electrons involved 
and K is their total momentum. 

For the ground state, single particle states are 
filled to the Fermi surface but the electron-electron 
interactions described above seem to be weak and varying 
slowly over the Fermi surface. In addition the energy 
involved in the normal-superconducting transition is very 
small. Cooper made a guess that the correlated states 
will involve single particle excitations in a small shell 
near the Fermi surface, figure 2.2. If we restrict our- 
selves to pair correlations, we shall want to know the 
interaction,which takes a pair of electrons from one state 
to another state in this shell. The total momentum of the 
pair is conserved and the phase space available for tran- 
sition from one state to another state of a given momen- 
CUM 12S a SELONGeLUunctiOn OL Ene sctoCal momentum, 7k 7gand 
isa maximum when K = 0. It is found that the effective 
interaction between electrons of singlet spin is stronger 
than that of triplet spin due to the exchange terms in 
the electron-electron matrix element. If all the pair 
correlations have the same total momentum, the maximum 
correlation of the entire wave function can be obtained. 
Fors dynamical and statistical reasons there 1s a strong 


preference for momentum zero, singlet spin correlations 
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Figure 2.2. Two intersecting shells each of radius 

k,+ 6 and thickness 6. The shaded area is a cross- 
section space phase available for scattering, with 
conservation of total momentum K, of a pair of electrons 
with individual momenta restricted to the shell 


kp =O. <.k S kt 6. 
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and for statistical reasons alone there is a strong 
preference for all the correlations to have the same 
momentum. For these preferences the wave function 4 
to represent electrons with the lowest possible poten- 


tial energy is made up of wave functions of the form 
UCP +b (=P +) 2.14 


where p(Pt+t) describes an electron with momentum P and 
spin up and ~(-P+¥) an electron with momentum -P_ and 
spin down. @’ is then represented as a total sum of 
these functions 

O (ry Dn) -= ) a; W(P,+)p(-By) REE 


where 


gives the probability of finding the electrons at any 


instant with individual momenta Pie 


(rj ,r,) describes what is called a "Cooper pair". The 


total energy of the two electrons is 


W = potential energy + kinetic energy 
2 
+ AE ae} 
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The states available for the two electrons are those 


with P> Pie EE ®(r)/r5) is normalised, 
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and the kinetic energy > 2P2/2m = 2€p where suffix F 


refers to the Fermi surface. 


(1) The Ground State and the Energy Gap 

mnie electron-electron interaction 1s attractive 
wnen theirvenergies e(k), e(k-q)Pare nearly equal’ and 
the matrix element is negative. The interaction is 
repulsive when their energy difference is less than 
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and the matrix element is of the form 
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The BCS criterion for superconductivity is that 
the phonon interaction dominate over the Coulomb inter- 


action and define the effective matrix element as 
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where K OFF is the effective momentum exchange. 

The electron pairing idea is extended to a large 
number of electrons and in the ground state all the 
Cooper pairs should have momentum zero. The ground 
superconducting state will be formed of those paired 
states which have the largest number of transitions 
among them. For pairing between electrons with equal 


and opposite spin, the reduced Hamiltonian is 
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where T is the kinetic energy term, Het is the inter- 


action energy term, Vick! is the interaction matrix ele- 
ment, a, and a, are creation and annihilation operators. 


The BCS ground state is given by 


= leet 
Ve ee z (u, + vp a, ai, ) |0> Pea 
where |0> is the vacuum state with no electrons. u, and 
Vv, are variational coefficients satisfying the normaliza- 


tion condition 


; 2e4 
inethis state there are no unpaired electrons. Vv, is 


the probability ef finding a pair in the pair state 
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(k+,-k¥) and uy is the probability of this pair state 


being empty. [25> becomes the normal Fermi state when 
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Ws is a maximum when 


(2.27) and (2.28) give the results 
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AL. turns out to be the energy gap parameter for the 


kth electron. Introducing another notation 
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Replacing the summation by an 


of states N(E) (2.30) gives 
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integral over the density 


where N(0) is the density of states at the Fermi level. 


Integrating (2.32) gives 
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Figure 2.3 shows’) the behaviour of equation 2.35. 
Experimental results on tunneling by Giaever et al 
(1962), Adler and Rogers (1963), Rowell et al (1963) 
are in good agreement between measurements on Al, ee 


Sn and the BCS density of states. 


(iii) Critical Temperature 

The energy gap is not a constant but decreases 
as ther temperature rises. An) electrongin stateg(P?t) 
without 7a partner in State (—-PY) prevents=the pair state 
(Ei efron being available yton the Cooper pairs and 
the pair interaction energy is decreased because the 
number of scattering events in which they can participate 
is decreased; and this implies a decrease in the energy 
gap. As the temperature rises, the number of free like 
broken up Cooper pairs i.e. quasi-particles, increases 
and the energy gap continues to fall until a temperature 
is reached at which the energy gap is zero. This tem- 
peratures Called the Critical temperature To: 
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d, Wipe Uber eltese spy tes bes Pha ay 2.36 


where 


bi 
T 
S 
7h 
oh 
mM 
nw 
i 
I 
= 
© 
ke 
5 
ke 
I 
im) 


yO9,-,L4 m0 etnonsauesiem disstdod | 
— 20 aes os edd bas a8 


Su hae 
ctw tion fsoista® (454) 


asussiosb gud tnsetesoo, 6 ton ei gs6pP yeusas oT 


(+9) .o3s2e at anortoslo aA .2028i% Sivtsiogmes eit 25 ; 


etase ting eft etasvetq (+97) stese oi tectisq 5 tuodtiw 
brs exisd YegooD sdt os efdsiisvs puted moxd (a= 4) 
git sausood beesoiosb ef ypisime soistosxesnk tisq ens 
eteqividxég aso yedt doidw of stneve paisedssoe xo tedmun 


yptene sds oi saaeior s eotlami etn bas vbsessxab my 


g#il sett to t6dmpn oid ,abart sus request arts BA .qse 


egasstoni .2eloisusq~iasup Lore ot bee 18qo0? qu nsaoid | 


suutexeqmes 5 Litay Lisi od aouattnoo gsp yptene eft bas 
-ste3 eifT .o19N 2i GBP ypT9a9 oad rio ite +6 ‘betosex ak 
“aT exusstsqned Sspiaito dd batino ei sxudeteq 

asvite (96.8) #0 cotssutsva 7 


af .S (egos - 00 PEG cr gill av Rf 2 A 


oterlw 


Prac tne 


. (S3eny 


al oe 7 | 


19 


Puc ulgem 2 


Density or svaves in the BUS theory: 


G. $ squat 
.¥iosdt oa ens mt ga date 10 fatten 


20 


e L. temp 
As Sen yin” Pp amass, 
e =r alk 


vsethetrerml function.) Substituting for h, we get 


Ke 2 
i a a 2B (7 E, /kT ) 
(e rend) 
putting 
A. = Ayes = A 
Vik! = V 


and using the integral instead of the summation, 


Aw Tes Ay? 
oomtanh. |) 
eee sols te siamese AS caer Jp de I 236. 
NCO)Vi © 2 Dee ¥ 5 
0 (estes) 


(2530) is equivalent to (2.32) at T = 0°K and determines 
iO) mentale eee Nee wO, ond for -s T es, (2°36) (has ino 


solution but the metal is normal with no energy gap 


Aw tanh Cpa! 
Cc 


l re 
TOyV = | scarica de EA accel 6) CA AR Se 
0 
For weak coupling, 
men ASIAOUNE 2.38 


j Se ane ey aie ae 
Cc Cc 


PLO (2.54) tallamtzicco) ace s=1Oln, 


2AN(0) = 3.0 kT. : 232 


) 
» Foe 
a D 


asnimssso5 bas 190 = 236. se. s) os Peete ns ea an 

on abit (QE-8) 4 gs < uot bas 0. + 4 a oe “a | 
(ia oe 4 

aoe! werent on ddiw ee pireses stig aud aolsules | . ~ 


) hei, ‘Sy dail aa 


| (40 = te +: cee.” aie 7 
or eee ee a 


: 
, : 
; ' oa an - 
7 a ae 2p | 7 Ris eel : | | | 
f 


y - : : 
i] _ re aay, | " 4 Po a 
hos. 


a %e Dans Pn es Ava A | | | ; 
vice rele me an A TG ree ee ee | 


Near the critical temperature To ; 


L 
NOB) = EPR (eh = =)? : 2.40 
(e 


The energy gap as a function of temperature is shown in 


figure 2.4. 


(iv) Cay mines eCriticalarield 


The critical field satisfies the relation 
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where G (0) and G, (0) are the Gibb free energy densities 
of the normal and superconducting phases. At T = 0°K, 
COU G, (0) is equal to the difference between the 
internal energy densities in the normal and superconduc- 


ting states. According to the BCS theory this is the 
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total binding energy of the Cooper pairs. Equation (2.41) 


then reads 
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Figure 2.4 
Temperature dependence of the energy gap in a 


superconductor. 
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(b) oize ELTfect) om the Critical Field 

When one of the dimensions of a superconductor is 
Otetne order of ethe penetration depth, thercritical mag- 
nNetveerieldsismmuchshigher than that of the bulk speci— 
men. 

In an external field He , a Superconductor acquires 


a magnetization M(H.) and becomes normal when 
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e 
fo mai.) au = 
e e 
0 


Equation (2.43) represents the area under the curve of 
figure 2.5 in the case of a spherical superconductor, 
Gynton 9(19162)) = 

Dhevsusceptibimlity y determines the initial slope 
of the magnetization curve. A lower x implies the curve 
HasmcOlGOmlLOmasnhigneryCchitiical field Hy for the small 


sample. If the curve remains linear up to Es 
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Because of the addition terms in the free energy equation 
(2.41), the penetration depth increases and the critical 
field becomes much higher. For a thickness of 2a, 
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CHAPTER 3 


MODELS OF TUNNELING 


Giaever (1960a,b) demonstrated experimentally 
that a current can pass between two metals separated 
by a thin insulating layer. This work was extended 
by Nicol et al (1960), Giaever and Fisher (1961), 
Giaever and Megerle (1961) for the case where one or 
both of the metals are superconducting. 

The foundation for the theory of tunneling was 
laid by Bardeen (1961) and further treatment has been 
made by Cohen et al (1962) and Prange (1963). Adler 


(1963) has given a good account of these theories. 


A. Simple Model of Tunneling 


If two metals are separated by a thin insulating 
layer, figure 3.1, electrons can tunnel through from 
one metal to the other. At low voltages, the tunnel 
current is ohmic. In this picture the transmission 
coefficient of an electron to go through the potential 
barrier varies as exp{- at (mo) 2} where m is the electron 


mass, nN Planck's constant, t~the-imsulating layer thick-— 


ness and » the work function. 


fealiieecacees ee unan (dys0808). xavesiD a - 
petaxeqee elsiom owt aaawsod aésq 169 sass 8 tens 
bebnesxs 2sw Az0W aide stoysl entsetuent pais 6 a 
, (haed) age bas xovesid ,(90@l) Is 39 Loot yd 
_%0 sto siodw peso aid 103, (L3@L) olxepsM bas revasia 
| -pnidsubnosysque o7a alsiom add to dtod 
asw poifsnaut to vioerl+ edt 10% noitsbavot sit 

asad asd Jnentsots startles ; Bas (Laer) nsebisa yd bist 
<SIbA .(€82L) sptisxa bas (S8@L) Ls Je asred yd ebem 
“eetioerts sesdji to tauovds boop 5 rey a esc (Ed@L) 


pofitalueni nkdd 6 Yd betsisqse sie elstem ows 21 
mort rpvorit [ennut nso enotiosfe \1.£ suupi? ,teyel 


Lonawt oft yaapediov wol tA .xerlto eft o¢ Isdom sao 


noleeimenisi1s ods sirusoiq eirbhh at .nimio at a 
{sittnedog od3 dpuorsdt op oF aoxtoole fs to dnoisittees 
soxjosis els ai m oxeriw £ com) 28. -}qxe 2s eorisv 1itisd 
-duidt xasysl pattslwent ent + (Ste seMOo @' Aone ld a ae. 
ote xxow ofd > bre sil 


q | 


a 


et 


eaten 
MF iigh 
T INSULATING LayeR 


QQ 


<—— -9- -——> 


1 bea UNOCCUPIED STA ee 
Ve (ZA OccuP\eD STATES 
E 


ro) 


V 
netinale ef ah 
Tunneling between two metals seperated by a Chain 


ineul@ting layer. 


28 


The spsobability that an electron will tunnel) from 
the sleitr metal to the right metal is proportional to 
the eyelet Mage fie that there is initially an electron 
coe riewlert metalinestate kk. This electron can tunnel 
POptenGemoht metal sinto state K'bonly arestate kis 
inittvally empty. —According to Fermi's rule, the tran— 


Sitzou probability is Given by 


ze 


P = 
kik, 7 


2 
|| fin =f NZ Byaal 
where M is a matrix element for the transition, ft. and 
ft. are the probabilities that the initial and final 
states are occupied and Ny TSoethemeinaledensityaotmectatess 
Tie .opam tial current. 1 mike eoLODOntLLOna LaLco 


alse 


2 
. Le N,N,£,(1 - f£,) |M| 
=i/=f 


where Ny and N, are the density of states in the left 
and in the right metals respectively. We assume that 
the tunneling electrons impinge normally on the boundary 
and that there is a one-to-one correspondence between 
the specification of states in ihe lert metal and, right 
metal. 


The net current across the barrier when voltage V 


is applied is 
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ie | IM|*{N, (E) £(B) (1- £(Etev)N. (Etev) 


= 00 


7 N, (EteV) f (Etev) (1 - £(E))N, (E) }dE Sie 


where E is the energy measured from the Fermi level and 
f is the Fermi function. 


When eV<< E,, , where E 


: is the Fermi energy, M 


1 


is taken to be independent of eV. 
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It is further assumed that near the Fermi level the 
density of states is constant. Equation 3.2 can be 


written as 


(oe) 


a ane jm? | N, (E)N. (EteV) {£(E) - £(E+ev) }dE 
o) eV 
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= “fs im|? | N,(E)N,(Etev)dE . 348 
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where in the normal state Q, (E) = Q. (E) = 1 and also 
eV is very small. Now the current can be expanded about 


the Fermi level. 
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Figure 3.2a shows the densities of states at zero and 


non-zero temperatures and the I-V characteristics. 


(b) Normal-Superconducting case, I 
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The corresponding density of states and I-V characteris=- 


tics are shown in figure 3.2b. 


(Ge) superconducting—super conducting case, ia 


In this case Q, (E) and Q. (Etev) are no longer simple. 


The current Tos is given by 


0 
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eV = 2h eV - 2A 
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where K and E are complete elliptic integrals of the 
first kind. Adler (1963) has carried out a rigorous 
derivation and calculation of (379). The density of 


states and I-V characteristics are shown in figure 3.2c. 
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At 0°K Tos is zero until the applied voltage is equal 
le) Aoaat A. at which voltage it will increase sharply 
approaching Ion at higher voltages. At non-zero tem- 
perature, a small current will flow at small applied 
voltages due to thermally excited electrons above the 
energy gaps. This current will increase up to V= Ao- Ay: 
(A, > A,) and gat Vv > A,- Ay it will decrease because the 
same number of electrons will be available for tunneling 
as in the V-< A,-A 


2 al 


+ Ay and at that voltage the current will rise 


range. This decrease will continue 


eiphes dL We SS A. 


sharply. 


B. the Modified Tunneling Model 
The effective Hamiltonian is given by 


H = H, + Ho + 4H, Sigeh®) 


where Hy and H. are the Hamiltonians for the left metal 


and the right metal respectively and Hn is the electron 


transfer Hamiltonian given by 


a (27 a) 
By ge hectare acts je sing 


a kk's 


where s refers to spin. 


Schrieffer (1963) considered H,, as a small pertur- 
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when the electrons potential decrease by V in moving 
TOM uatO 25 |o4> and |B5> are eigenstate vectors of 
Hy and H respectively. Ei and Ep are energy diffe- 
rences between |0,> and Ja >. |0.> and | 85> respec— 
tively. 

At 0°K electrons cannot tunnel from 2 to 1 because 


of the conservation of energy. The current I(V) is pro- 


portional to the rate of transfer of electrons: 
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(+) and (-) refer to "electrons" and"holes" respectively. 
Using Green's function, Schrieffer reduced (3.14) and 


(3215) “to 
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Differentiating (3.13) with respect to Vv, for the normal- 


superconductor case, 
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Gis called) the \dynamical conductance. From (2.35), 
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The process thus described is illustrated in 

Pi gusess 2 -— = tnewonayan=eclectrongin kt) below=the 

Fermi level in the normal tunnels through the barrier 

to state k'+ above the Fermi level in the superconductor. 

A hole is left behind resulting in an excitation energy 


e = |e 


a for metalgl. Due to a quasiparticle in k*4, 


«| 
there is an excitation energy hs Eyes 
in metal 2. This process will occur only if the pair 

State (kit; ee is ifitially empty?) Energysis conserved 


if je, | + E,, = V- In figure (3.3b) the particle goes 


in state k" below the Fermi level in) the superconductor. 


if A. is constant, Gy. tame ae ue 
In figure (3.3c) where both metals are superconduc- 

Ging paancelectronvin Kt Gan> tunnel ito vei ther kotor 7k 4, 

or an electron in kt can tunnel to the same final states 


The current this time begins at V = Ay + A. 


Guam the Capacitor Model 


In an attempt to explain the observed zero-bias 
resistance anomalies, Giaever and Zeller (1968) intro- 


duced a capacitor model. Mezei (1971) has tried to 
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explain this model in a microscopic picture. Townsend 
et al (1971) have followed Giaever-Zeller argument 
except that they discussed the dependence of conduc- 
tivity on zero bias and temperature at both zero and 
non-zero temperatures. 

The potential energy of a sphere of radius r and 


carrying charge e is 


2a 2rerey 456 i 
re) 
where Q = e and the capacitance C = Amer. Li contact 


is made between the two sides of a parallel plate capa- 
Citor, figure (3.4), which has) capacitor plates made of 
metals with different work functions, electrons are 
exchanged until the Fermi levels are aligned. If one 
Op the plates 1s oL va small tinite size, exact align- 
ment is not possible; alignment occurs in discrete 
voltage steps AV = e/C. At equilibrium the number of 
electrons on the particle is closest to alignment, the 


difference in voltage being |V < 6/2C. At low tempera- 
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tures the number of electrons is fixed. The two films 

and any one particle can exchange electrons, figure 

(32 5a) through tunneling untilvequilibrium is attained: 
When the voltage is applied, electrons can flow 


from one side of the junction to the other side by any 


of the following processes: 
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Figure 3.4 
Particles embedded in oxide seperating two metal 


films represented as two Capeci lors iio ceies. 
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it Direct tunneling through the metal oxide as in a 
usual metal-metal oxide-metal junction. This gives a 
constant voltage temperature independent conductivity. 
1 Ttunne lings trom one’ film to a particle lecalising 
these leéctronstheme, for a short time andythen tunneling 
Gut eo the Other faim. Whatetollows will be above this 
process. ie 
III. Tunneling as in the second process but not loca- 
lising the electron. 

Process II requires an activation energy which 
in turn is responsible for zero bias resistance peak. 


The activation energy E is equal to the energy difference 


Detweenesthe. internal, and, finabistates. 
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The level spacing in the particle is assumed much less 
than E and therefore the spacing effects are negligible. 
E is purely a classical Coulomb energy. 

In the pictorial representation of figure (3.4), 
the total capacitance is 
ae Peau 
CR Ge Cy, 


C' = 


Gh M.,-oxide 1s much thicker than M,-oxide, Cr >? Cp so 


that the applied voltage appears between the particle 
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and the right film. “When alpanticlée munnels from left 
to the particle, charge flow is opposed by a polarising 


charge 


due to the left capacitor. The net charge is 


dl ec, ° eCp 
Cr t+CR C.+CR 


Battery energy is 
ec, V ec, V 
qv = = . 


Cl tCp e 


For, currenteto, Llow the applied voltage must be greater 


than the activation energy, i.e. 


evc 2 
L Ane 2 Lael 
Cae. tC ele Wb ec oD 
C S 
V2 aq (ga # Vp) . 
L 
With the assumption C=C, >> Ce pacDeRCONCItLOn £Lonecurrent 
flow becomes 
lv} >ortV fon ee a9 
XE D 
\vjo> sev for eae ay 
2C D 


Therefore for a single particle, 


ai ypisne yrestsa a 


- 
io ae 
Vv Vewe!. . = a= 
. = =< - 5 = = Vp > 
| a 7h 
vu - 7 
: 


tetse1p ed teum spstlov boifggs sid wolt ot JneumuD SOT 
te .9.i ,ypasds notisvisos sat asds 


On (V) Opaest - 55+ Vp <V<aatVy 
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constant otherwise . 


Figure (3.5b) illustrates equation (3.19). We assume 
that for a set of particles of given capacitance C, 
Vp is uniformly distributed in the range -e/2C <V<e/2C. 


The average conductivity for this set of particles is 


e/2C 
ome ay e/2¢ 
aca Miies eh See re ey ae o (Vv) av 3820 
by ee72c 7 eS D D : 
dv, -e/2C 
-e/2C 
Define 
ee 
V4 = 5a + Vp roost 
(3.20) becomes 
e/C 
C 
SOV a= oa | Op (V) qv Smee 
0 
and On (V) = 0 ie ve Se ON 
e 
or ie < Maes 
On (V) =o otherwise . 


1 is chosen arbitrarily as any constant. But the limits 


for integration are zero and e/C. 
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For Ueca.Vi , alae Vie = ; 
e/C V 
| Op (VV) dv. = | qv, = V 
0 0 
. e 
gia Vane ral 
e/C e/C 
= ka 
0 0 
For Os Vie, sie Vo- = 
e/C -e/C 
| On G4.) C= f av. = -V 
0 V-e/C 
& 
hee V<- ral 
e/C e/C 
= ae 
| Op (V,) qv | av_ = a: 
0 0 


Putting these results in (3.22) gives 


O|Q 


<o (V)> |v | for Ivl<=s o.26a 


o1AQ 


<a (V)> -& = 1-= constant for |V|>& 3.23b 


The variation of <o(V)>" with voltage is “shown an figure 


(356). 
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Figure 3.6 
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YoOrravdistri bution OLeparer cles seuch that. £(C) 


is the number of particles with capacitance C, 
o (V) =| (elere (Sy) ee VAS A 3.24 
0 : 


The factor C comes in since o(V) depends on the capaci- 


~ 


tance which in turn depends on the area of the capacitor 
plates sapUsing: (oF 23hmwand takingsposselVverV; 


e/V 00 
face*ec) =F | CaGEtC)a a: S504) 
0) ay) 


o(V) = 


OI1< 


The number of particles f£(C) is determined by the dis- 
PEUDUELOneOLspaltt Cles sr lice CONnGUCLi Vit yeg iy ene by 
equaclon so .25e sel Lustrated sin figures .600))) 14) 
increases linearly for small applied bias voltages and 
the main contribution comes from the large particles 
Since for these c*£(C) is maximum. 

Townsend et al (1971) have discussed the case for 
non-zero temperatures. Their derivation gives the first 
term of equation (3.25) > temperature dependent. 

For particles of given capacitance C, the number 
Fi(C)mol excited particles is) proportional toe exp (-eV ,/kT) 
and the number of particles of capacitance C with activa- 
tion energy between We and Votdv, is . £(C)dv.. Thus 


eV 
Cc 


F(C) = £(C)dv, exp (- ae . 
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The contribution to the conductivity is 
c? =-eVA/ kL 
=e E(C).d¢G e Cc 
- ne 2 -eV/kT 
o(V,/T) = | = EXC) aC | e av 
0 0 
fee) 
x <. | c*F(c) act 3.26 
ri tc 
ac Vv = 0 , 
c 


The second term which appears in (3.25) is left out in 
Use 20)". 


e/C. 


C(O TD) er 3 


When the particles become superconducting, the 


It becomes important at voltages greater than 


46 


activation energy increases by half the energy gap; ev, 


becomes ev, teh. Equations (3.19) las tor pe rtepiaced 


by the normal metal-superconductor voltage dependence. 
These changes will cause (dI/dV) V9 


well be“reflected ian “‘the™ conductivity as the critical 


temperature is reached. 
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CHAPTER 4 


EXPERIMENTAL PROCEDURES 


ee Sample Preparation 


A tunnel junction is simply a metal-insulator- 


metal sandwich , henceforth referred to as M,-I-M,, 


deposited on a suitable substrate such as a glass 


microscope slide. Such junctions with particles 
Gepocttedwinethesincscularor,;sreterredslOeds M,-M5-M, 


were fabricated. 

A glass microscope slide was cleaned with a 
detergent under running tap water, dried with optical 
Wiper tissue and then passed over a hot flame to remove 
impurities by ionization bombardment. 

In the case where M, and Mo were tin, indium metal 
Ben cncts were made on the glass microscope slide before 
mounting in the slide holder. If contacts were made 
SpteuecvapOkaclOnwoOr, tin, ALONG Lit sti lm wacm~OUnNd CO 
break due to the low meLeinG DOM Oi et li. 

Evaporation techniques in the preparation of thin 
films and particles have been discussed by Holland 
Glo oy mand, CoOpram | 0969). sALtoreatlcalningeasvacuum 
OF insu Tome torr in the evaporating—bell—jar system, 


M, was deposited and the surface oxidised as described 
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in the following section. M, particles of different 
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sizes controlled by the quartz crystal thickness monitor 
were deposited and then the surface was allowed to 
oxidise slightly. Finally M, was deposited to complete 


the junction. The stages are shown in figure 4.1. 


f)-Oxidation 

Different methods of oxidation were tried. 
(i) Aluminium 

Aluminium oxidises easily even at room temperature. 
Aluminium oxide of desired resistance was formed by glow 
discharge under 0.1 torr of oxygen and a total current 
of 1 mA. By allowing longer times of current flow, 
oxides of higher resistances could be obtained. With 
the aluminium film 4.5 cm from the aluminium electrode, 
8 minutes duration gave a tunnel resistance of about 
200 2. Miles and Smith (1963) have given a good account 
of this method. 

Thermal oxidation in air under humid conditions 
and in dry oxygen at one atmosphere were also tried. 
In the latter case the glass microscope slide was held 
over a heated semi-cylindrical copper sheet. The tem- 
perature was recorded by a thermocouple. Again by 
fixing the heating current and hence the temperature, 
the duration time was the major factor determining 


the resistance of the oxide. 
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Cid) maar 

Tin is more difficult to oxidise. Good barriers 
of tin oxide were obtained by evaporating tin at a rate 
of about 4 A per second and oxidising the film surface 
by thermal oxidation for two to four hours. 

Tin forms continuous films at relatively thicker 
films and the film structure depends on a number of 
things such as the nature and temperature of the subs- 
trate, the temperature of the evaporating source, rate 
of evaporation and the residual gas pressure in the 
evaporating system, Chopra (1969). 

At high rates of evaporation, tin molecules impinge 
on the substrate in large numbers, and nucleation and 
coalscence of the grains proceed rapidly forming large 
islands and then a milky continuous film. Such a film 
is difficult to oxidise, Walmsley (1965). But at low 
evaporation rates, no milky film is formed and the film 


is easier to oxidise. 


(iii) Lead 

Lead is easily oxidised in an atmosphere of oxygen. 
With aluminium as a base electrode, the tunnel resistance 
is much larger but it does not increase with time as in 
the case of Aluminium-Aluminium oxide-Tin. This obser- 


vation agrees with Handy's (1962) results. 
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Oxides of lead were formed by exposing lead 
films to one atmosphere of dry oxygen for 5-20 minutes. 
In the case of lead particles, the resistance of the 
tunnel junction increased rapidly with increasing 
thickness of the lead particles. It is suggested that 
with the low melting point of lead, for small particles 
a number of them have the oxide over them breaking down 
on evaporating the final electrode, aluminium,which has 
a much higher boiling point. As a result the resistance 
is reduced for very small particles and one gets a 
situation with Al-I-Al, Al-I-Pb-I-Al and Al-I(Pb)Al1 
altogether forming the tunnel junction. This suggestion 
is supported by the results of junctions with small par- 


ticles of lead. 


Ti) DepositionNofeParticles 

In the film formation process, after nucleation, 
the grains coalesce into small particles which then 
amalgamate to form a continuous film. The idea is to 
get amounts of evaporated metal M3 without reaching the 
continuous film formation stage. 

After a constant rate of evaporation of metal M3 
was attained, particles were deposited onto the oxide. 
The rate was about 10 A per second in the case of tin 


particles and about 20 A per second in the case of lead. 
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The amount evaporated was controlled using a quartz 
crystal thickness monitor described in section B. 

Ines size onde distuasoutYonwote partie cs Ot eleadmwas 
observed with the electron microscope, plates 5.1 and 
Serer: 

The surfacesof the particles were oxidised to 
form achat oxide. © ini the case of tinsyathermeal oxi — 
dation was for about two hours while for lead a few 
minutes in dry oxygen gave a reasonable oxide thickness. 

Particles for observing in the electron microscope 
were deposited on an oxidised Al base held on the elec- 
tron microscope grid by surface tension, simultaneously 
with the junctions. The electron micrographs were 


analysed as described in section B of Chapter 5. 


B. Thickness Measurement 


The thickness of an evaporated film was monitored 
during evaporation by a quartz crystal thickness monitor. 
Two 5 MHz crystals mounted in the same water 
coGled Holder, Ligure (4.2c), One exposed tO, uie ecva-— 
porating sample and the other in thermal contact with 
the copper holder for heat compensation during evapora- 
tion were connected to the quartz crystal monitor as 
shown in figure (4.2b). 

This kind of thickness monitoring device utilizes 


the principle that the resonance frequency of a piezo- 
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Figure 4.2 


m Quercz Crystal Whickness Monitor 


(a) The Ssctliation of “a singe crystal; most. of the mass 
displacement occurs at the crystal surface. 

(o) The circuit diagram of the crystal monitor used in this 
work. The integrated circuit is a Motorola Me TG) 28 ee, 

(c) Crystal holder; very light springs are used to provide 
electrical Contacts with the crysvalse. 
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electric crystal depends on the mass of the material 
deposited on the surfaces of the crystal, Eschbach 
and Kruidhof (1965). The oscillation frequency is a 
decreasing function of the mass of the film deposited 
on the crystal and changes quite linearly with film 
thickness as long as the crystal is not overloaded.. 
The rate of deposition is found by observing the 
change in frequency Af, which is proportional to the 


mass added, Am 


Shy See 4.1 


where m is the mass of the quartz crystal of natural 
frequency fo: When the crystal is excited to thickness- 
shear oscillations, figure (4.2a), fy is related to 


the thickness, d, by the equation 


where C, is the transverse velocity equal to the velocity 


dh 
of sound in quartz (3340 m/sec). The frequency shift is 
given by 
or 26 
Af = - —sAd = - eee Ad 
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where PG is thesdeisity/oL avquartz ecrystal with area 
A and Am is the change in mass. If Am is the added 


mass of the evaporating material with density op, 


Me = 
hh 


where t is the thickness of the evaporated material 


Cie 
t=— 2 at . 4.4 
Z2Ee0 
fe) 
The crystal can be lightly damped at its periphery 
without disturbing its mode of oscillation appreciably. 
Warner and Stockbridge (1963) found that this descrip- 
tion of the crystal monitor is accurate within 1%. 
Hartman (1965) found the density of thin films to be 
less than that of the bulk sample. 

In our case, experimental arrangements are such 
that particles from the evaporating sample strike the 
crystal at an angle 6 with the vertical, so that there 
USeagcactor (OL cos 76 in the constant Bae cae For 
Stidial angles cos 6 “i. Since thes film density is some- 
what variable, the constant is best determined experi- 
mentally using the crystal measurements together with 
the Tolansky interferometer measurements. Vi Dean noi) 
has described these measurements in detail. The values 
used in this work are shown in table (4.1). In all 


thickness measurements, t will refer to average thickness. 
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Tables 431 
yds 
Material o (bulk 2 o- = ppm/KR 
gm/cm 1°q 
Al Ze70 220 
Sn ae 75 440 
Pb 1ik 34 850 


Cs. iWroeduction of Low Temperature 


The experimental arrangement for the production 
of low temperature is shown in figure (4.3). In the 
figure, (P) is the pumping chamber provided with a 
needle valve (a) to allow in the chamber the desired 
amount of liquid nitrogen or liquid helium. A stainless 
steel tube (b) provided with thermal shield (h) connects 
(P) to a rotary pump. By pumping over liquid helium, a 
temperature of about 1°K could be reached. A regulating 
pressure valve to control the pressure above the helium 
allows a variation of temperature between 1°K and 4.2°K. 
Elevated temperatures up to 70°K could be maintained 
using an electronic controller which senses the resis- 
tance of a 220 2 speer carbon resistor (d). (s) is the 
sample chamber which contains the sample holder (e) and 
the germanium thermometer (g). The vacuum chamber (V) 
is connected to a diffusion pump and a pressure of 


6 


1x10 ° torr is needed to shield the sample from the 


liquid helium bath (t) which is surrounded by liquid 
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Figure 4.3 
arrangement for the OroductLlon 
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nitrogen bath (N). The cryostat was designed to fit 
in a bore of a 20 KG superconducting magnet (M). 
Indium 0 rings (0) plus apiezon grease were used for 
vacuum sealing. The grease makes it possible to 
separate the cans (c,) and (c.,) without damaging the 
indium 0 rings, thus making it possible to use the 
seals four or five times before repairing them. 

For temperature measurements the germanium thermo- 
meter #904 was calibrated by Dr. Rogers between 1°K 
and 60°K and by Dr. Woods between 4.2°K and 12°K. The 


bridge used was designed by Dr. Rogers. 


Dis Conductance Measurements 


The dynamical conductance o = di/dV and its deri- 
vative do/dV = d*i/av are measured as functions of the 
bias voltage V. The conductance bridge used for this 
purpose was designed by Dr. Rogers and has been described 
by shim, @Rogers (1970). TA Lock iurcamphities Iseused tas 


the detector and an X-Y plotter for recording the results. 
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CHAPTER 5 


EXPERIMENTAL RESULTS AND DISCUSSION 


Aw TinsParticles 
(i) Sn-Sn’-sn Junctions 

A conductance dip at zero bias is observed even 
at room temperature but more pronounced at 80°K and 
below. JOr the three junctionseshown in figure (5.1 )y, 
the one with average particle thickness of 50 A has 
the greatest dip. 

In figure (5.2) the dip is very large. At low 
voltages the conductance-voltage dependence is linear 
but rolls off at higher voltages. The general shape 
of the curves is similar to that predicted by the 
capacitor model, figure (3.6). 

At low temperatures when Sn is superconducting, 
figure (543), a peak at JT 7l5SemV and a step at abcut 
2.5 mV are observed. The peak decreases as the particl 
size increases implying that direct tunneling between 
the filmsjjis_not—important—forlarge«saze-partieless 
The position of this peak determines Che energy gap of 


Sn where V = 2A Tunneling between the films becomes 


Siw 
important as the particle size becomes ,smaller.,...This 
is so since the activation energy for smaller particles 


is higher. 
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ARBITRARY UNITS 


O 50 100 150 200 


V(mvV) 


Figure 5.1. Voltage dependence of the conductance 


o of three different junctions of Sn-Sn’-sn ate BOCK. 
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Figure 5.2. Voltage dependence of the conductance o of 


three different junctions of sn-sn?-sn at 4.20K. 
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2.1°K 


t= 20 A (a) 
50 A (b) 
65 A (c) 


Vv (mV) 
Figure 5.3. The conductance yersus yoltage for three 


different junctions of Sn-Sn’-sn at 2.1°K. 
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The step occurs at about 2.5 mV and this is 
equivalent t : : imi 
q fe) 2Aciim ae a particles , samilar to two 
junctions in series. In section B we shall see that 


the second derivative of i gives a more precise point 


for this step. 


(ary) Al-sn*-al Junctions 

The voltage dependence of the conductance for 
Al-sn'-al junctions shown in figure (5.4) is similar 
to the Sn-Sn?-sn junctions at 80°K. At liquid helium 
temperatures, figure 5.5, the zero bias conductance 
anomaly is also similar except that the rolling off 
occurs in a smaller voltage range. For small bias 
voltages the conductance is linear rolling off at higher 
voltages depending on the average particle size. For 
large particles, for example curve (e) (t = 130 Aygeit 
levels off almost sharply. The conductance dip at zero 
bias, o(0), varies almost linearly for small particles. 
This is shown on the bottom right hand corner of figure 
(5.5). The half-width becomes narrower as the particle 
size increases. The shape of the curves agrees with that 
Or Ligue, (3.0) % 

At lower temperatures, shown in figure (5.6), when 
aluminium is superconducting, the energy gap cannot be 
noticed implying no direct tunneling between the films. 


‘The step at about 1.4 mV can hardly be distinguished. 


€3 


wie witag 5 na ale ole 
tent ese Ilsde sw a noitose al .aeixse ait enotsonut 
tniog seissiq siom 5 sevip i to evistsviwsb baovee eft 


.qste eift 10% 
‘ 
enottonut Linea (BR) 
tot songstoubnos sit to sonshusgqeb spasiov aT : 
aplimte af (b.2) edapha nb awode enolzonut ta~*n@-1A 
mu ifert biupil- 4A .1°08 $8 anotyonut n2~"92-c8 sits 03 
SonkdGeerios ebid ores eds ,c.¢ sinpLt , 29tu35 rEqmes 
4to pat{fox, sit tedd sqeaxs sgl Aachaa ol at yismons 
estd Iisnte 0% -Spnsz epasiov a1sileme 6 ak 214090 
rsipid ip 220 priflGs ~zeeaif eat go rte pron ens aspstfov 
40% .oxts ofottrsq Spsxeve sdt ao prbbitegeb sapasiov 
ti ,(A O€L = +) (9) svasp sigqmsxs x03 aalsitisg spisl 
otes +6 gibh sonstoubnds odT .ylgusde teomis ito eafsvel 
ealoisisg ifisme so? yineenis jaonis asitey ,(0)o ,esid 
siupit to sx9nto> Basi Jnpit medsod sft no Awode et allt, 
‘elotsisg eit e6 toworrsa samoved asbiw-tled off .(€.2) 
ject didiw aeoxps esviw9> eit Xo oqeda sclT _+S9esetont este. 
«(3.€) suvptt tom 
aadw is bas ceils nt a si iiiep ated towol 3A ; 


4 


Te iW a vo af sang! \ ats! os at; 
ee on =a ia *) a 


64 


ARBITRARY 
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Figure 5.4. Voltage dependence of the conductance o 
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ARBITRARY UNITS 


118° K 


25 A (a) 
50 A (b) 
65 A (c) 
908 (d) 
IZ0B (e) 


V (mV) 


Figures5.6. 4 Conductance-voltage characteristics of 


five different junctions of incisal atm laws tee 
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The voltage at which the conductance saturates depends 


on the particle size as can be seen in curves (e) and 


(cd) Withivervy thick particles showne1 pet iqu coms 7, 


two junctions in series behaviour becomes prominent. 


The first peak at 0.8 mV is equivalent to (Any a IN 


Sn? 


and the next one at 1.27 mV is equivalent to 2A and 


Cnew chide Olewatel9 MV COLresponas™=to (2A, 


it 


ae 2A..)- 


one would suspect this last peak to reach a maximum 


when the 


Pabticies form,a, continuous, fiim-: 


Figure (5.8) shows the conductance-voltage charac- 


teristics for different applied magnetic fields. The 


curve at 


zero field is displaced when a field of rea- 


sonable magnitude is applied. The maximum field that could 


be applied due to the superconducting magnet was 20 KG 


ana wer couldmnot reach the critical field. 


The temperature dependence of the conductance is 


shown in 


paaleieverce (Gyo SONG Mme alici AkahigiclcWe sleke Gelcuielsvels muha sie 


above the critical temperature ue OL Sn and this coniimms 


equation 


(3226). When Sn particles become superconduc-— 


ting, the conductance decreases moving away from linearity. 


When the 


65A particles 


= 04. 15.2 0X0 be K LOG e 


130A particles 


= 4.05 +0.05°K Loi re 


3/20 Kee tOGEDULK eS mel vaOn eh Oloee ne 


linear curve is extrapolated, the o(0) intercept 


gives the background conductance of 27%, 4% and 14% for par- 


ticles of average thickness 65K, 90A and 130A respectively. 
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ARBITRARY 
UNITS 


V(mvV) 


Figure 5.8 


The magnetic field effect on the conductance- 


voltage characteristic. 


ARBITRARY UNITS 


(a) 


TOU DEM Or. Ocak 


conductance oO; 


Us 


The temperature dependence of the 


average particle thickness, t=65 A. 
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(b) Taek 
The temperature dependence of the conductance 


Figure 5.9(b). 


Oo; average particle thickness, t=130 A. 
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According to Blumberg and Seraphin (1963) the transi- 
tion temperature of metals such as In and Al, and Sn 
in some respects, whose a2 increases with decreasing 
thickness To is connected by the relation 

— 


A 
Sy SEP rd” bod Sil 
A t 2 


where A and B are constants and t is the thickness. 

Our tT. results at least show that T. for Sn particles 

is higher than the value for bulk Sn. The temperature 
dependence of the normalised conductance Ga/cs is shown 
in figure 5.10 for three different junctions. For the 

90 A curve i = 4.5 +.1°K and one would expect this 

curve to lie between the 65 A and 130 A curves. This 
does not happen because according to Giaever's suggestion 
there is a maximum size effect somewhere between 50 A 
and 100 A for Sn particles. Below 50 A and above 100 A 


particle sizes the effects fade out. 


B. Lead Particles 
(i) Particle Distribution 

Particles and thin films are usually prepared by 
the condensation of atoms from the vapour phase of the 
material. The way in which particles and thin films 


come into being can be ascribed to this growth process. 
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Nucleation is the birth stage. It is a problem of 
vapour-solid-phase transformation. Pocza (1967) has 
discussed a number of parameters influencing the 
structure formation. Such parameters are M, the nature 
of the material evaporated characterised by purity, 
chemical reaction, etc.; V, the vapour beam which 
would be atomic or molecular, characterised by tem- 
perature, pressure and evaporation rate; P, the 
residual gas pressure; S, the substrate characteristics; 
and W, the evaporation source characterised by geometry, 
temperature, etc. All these parameters interact with 
each other making it very difficult to survey them. 
Chopra (1969) has discussed in detail the formation 
stages of thin films. 

The impinging atom is attracted to the substrate 
surface by the instantaneous dipole and quadrupole 
moments of the surface atoms. It loses its normal 
velocity provided the kinetic energy is not too large. 
Thus the vapour atom is physically absorbed, it may not 
be thermally equilibrated, it may move over the surface 
by jumping from one potential well to the other implying 
that the absorbed atom has a finite "stay" time on the 
surface during which it may interact with other atoms. 
Alternatively the vapour atom could bounce off especially 


if its kinetic energy is very high. The probability 
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that an impinging atom will be incorporated into the 
substrate is called the "condensation" or "sticking" 


coefficient On 


where the T's and E's are the root mean square tem- 
peratures and kinetic energies of the incident (I), 
reflected (R) and substrate (S) respectively. 

For two different surfaces,say Al50, on an 
electron microscope grid and Al metal on a quartz 
crystal, which receive equal numbers of vapour atoms, 
the distributions of the condensed particles may 
be different or one surface might have more vapour 
atoms reflected than the other. We assume that the 
distribution of particles observed in the electron 
microscope is the same as that in the junctions pre- 
pared simultaneously. Plate (5.la) shows an electron 
micrograph of Pb particles on Al,0.. The average 
thickness of such particles is 180 A measured by the 
quartz crystal thickness monitor. Plate (5.1b) shows 
the diffraction pattern. There is a high probability 
that particle distribution may not be uniform even 
within the same junction. Plate (5.2b) shows this. 
In plate (5.2a) the large dark spots in the background 


are assumed to be due to the detergent used in the 
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Plate 5.1 


(a) Electron micrograph of lead particles 
of average thickness 180 oe 
(b) The diffraction patternnots particles 
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preparation of the base for the electron microscope 
Sparks 

Particle distribution functions have been 
derived by Chakraverty (1967) and Basset et al (1959). 
The functions are of statistical nature. The electron 
micrograph of plate (5.2a) was enlarged ten times and 
areas of the visible particles measured giving a his- 
togram shown in figure (5.11). The dotted curve of 
figure (5.11) was approximated to the curve of the 
function 
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2 distribution with 10 and 4 degrees of freedom. 


£or the Ww 
Equations (5.43) was substituted in equation (3.25) for 
£(C) and numerically worked out on the computer. 

Figure. (5.12) shews the curve plotted from the calcu- 
lations. By, adjustinggthe units, the measured conduc— 
tance at 8°K and the calculated conductance are compared 
ine froune (52413)% 

From the mean area measurements, if we assume the 
particles approximately spherical, the average diameter 
was 110 A compared to the average thickness of 60 A 
measured by the quartz crystal thickness monitor. ~The 


tTatio.ofetetal areéasgot particles to ltotra loaresa Grlelec— 


tron micrograph was 13%. 
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Figure 5. LI. 


particles of thickness 60 A. 
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Figure 5.12. Computer curve for o(V) versus VOLGage 


characteristic for Pb particles of average thickness 
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Figure 5.13. Measured (++ line) and calculated (solid line) 


conductance of a set of particles in the normal state, for 


particles of average thickness 60 ce 
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(ii) aAl-pp’-al gunctions 

in contaceawrtch Al,0.,, lead formed increasingly 
high resistance junctions with increasing lead thick- 
ness even after a few minutes exposure in dry oxygen. 
thas) made wut disfticult to getuqoods junctions. 

Figure (5.14) shows the voltage dependence of 
Chesconductancesor fivesartferent, junctions ateo.L ok. 
Particles as small as 20 A thick have no noticeable 
effect. The particle effect increases rapidly with 
increasing particle thickness reaching a maximum par- 
ticle size effect and then decreases with increasing 
particle thickness. The shape of the curves is similar 
to that predicted by the capacitor model. 

in) LEQuULe mo .LoOmtie tous Unec Ons =eSnOWwsve Ly dn— 
teresting differences at 1.18°K. In curve (a) the 
first peak is at 0.35 mV equivalent to 2An and there 


dh 
is a small peak at 1.45 mV equivalent to iy we 


Pb* 

In curve (b) the peak at Any is diminished implying 
that direct tunneling between Al films is unimportant 
with increasing particle thickness. There is an extra 
peak at 1.16 mV which does not appear in curve (a). 
This peak is equivalent to Abb = Anal and the third peak 
in the same position as peak two in curve (a). These 


last two peaks suggest that there is tunneling between 


Al and Pb particles without localising the electrons on 
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1.18° K 
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Figure.) «15... Conductance versus, voltage of four different 


JUNCELONS. On Al-Pb’-Al at, L.1eeks 
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the particles. | inecurves’ (cj) and (d) tunneling via 
particles is prominent. The second derivatives of 
such curves are shown in figure 5.16. .tn addition, 


one sees 2A, Table 5.1 summarises the results. 


fy 


Tabie Se CL .08c) 
Particle thickness 2Any, Particle energy gap 
A mV (20,5 + 2A5)) mV 
60 Zee 2.94 
80 2n02 2500 
130 2.66 3.06 


Some junctions had multiple peaks at temperatures 
when Al was superconducting. Figure (5.17a) shows 
conductance-voltage characteristic when Pb particles 
are normal characteristic of particle effect. When 
both Al and Pb are superconducting, figure (5.17b), two 
peaks appear at 0.27 mV and 0.41 mV instead of one at 
-35 mV, and another small peak appears at 1.4 mV. This 
is a result of some particles being in contact with Al 
film giving rise to what is known as the "proximity" 
effect. Vrba (1971) has observed such multiple peaks. 
De Gennes-Werthamer (1963) theory gives the properties 
of proximity sandwitches. When a superconductor is in 


contact with a normal metal, the normal may become 
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Figure 5.16. 9second derivative of the conductance 


showing energy gaps for three different junctions. 
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Figure 5.1/ 
(a) Conductance versus voltage with Pb particles normal 
(b) Multiple peaks of superconducting Al 
(c) Second deriyative of curye (b) 
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Superconducting while the transition temperature and 
energy gap of the superconductor are.lowered; if in 
contact with a superconductor their properties may 
change resulting in energy gap and transition tem- 
perature values somewhere between the two. 

The temperature dependence of the conductance 
was found to be linear when the Pb particles were 
normal in agreement with equation (3.26). When they 
became superconducting the conductance decreased 
moving away from linearity. Figure 5.18 shows the 
normalized conductance versus temperature. For the 
three different junctions the critical temperature is 
found to bé the same,7.4+ .05°K,compared to 7.2°K for 


the bulk sample. The average ratio 2A/kT , was 4.14. 


(ld) Phonon structure 

In Chapter 2 it was mentioned that Frohlich 
discussed the attraction between electrons as arising 
from electron-phonon interactions. Rowell and Kopf 
(1964) studied the phonon spectrum involved in the 
interaction. Solutions of the energy gap equation 
by Swihart (3963) ana) Culler at al (1962)) Snow 9 that 
phonon structure should occur at energies A+E, and 
A+E,, where EY is the Debye cut off energy at cut off 


frequency Wo and En is tthe Einstein peak. Morel and 
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Anderson carried on further calculations showing that 
for the Einstein peak, the density-of-states structures 


should be observed at AtE A+2E,, ; A+3E., Pr sneel 


E v 
Adler et al (1963) have observed phonon structures 


Of lead, in and indiun.."From cquatten (3.17b) 
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But due to quasiparticle damping and life time effects 


one includes an imaginary part in the energy gap. 
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Adler et al define the phonon strength factor S(V,T) as 
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According to the Schrieffer, Sealapino and Wilkins? (1963) 
model, Ay becomes very small and As increases at the 


phonon voltages Vo: 
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Figure 5.19 shows the phonon structures and the results 


are summarised in table 5. 


Tab le eo 2 


Observed Phonon Structure Voltages 


Particle size (A) ee 


There is a broadening which smears out the singu- 
larities. This broadening is predicted by Dickey and 
Paskin~ (1968)"in' their Calculations of phonon, spectrum 
changes in small particles. Vp depends on the particle 
Ssize.w Lt. is higher for ssmall pantrcles and decreases 
tending to the bulk value of 8.5 mV as the sizes particles 
inerease . 

Table 5.3 shows the particle results compared to 


the bulk values. 
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CHAPTER 6 
CONCLUSION 


The capacitor model has been used to explain 
results of tunneling via particles in the normal state. 
In the superconducting state tunneling via particles 
exhibits two junctions in series characteristics. 

There is a particle size dependence of the cri- 
tical temperature of tin particles. Bucket and Hisch 
(1952) found that the size dependence of the critical 
temperature of a thin film depends on the temperature 
at which the film was formed. All our junctions were 
prepared at room temperature. According to Bucket-Hisch 
observation, a higher percentage change in To would be 
obtained if the particles were deposited at a lower tem- 
perature. Due to the low magnetic field attainable, we 
could not keep the particles normal at low temperatures 
by applying a magnetic field. 

The energy gap results do not conclusively show a 
particle size dependence of the energy gap. In the case 
of lead particles 2hap show,.a,2.2%, variation about, 2.65 mv. 
TiereslSmasbDigGguncer taut yin Detios ab LO 2A/kT , since To 
could not be determined precisely. 

The phonon. structure in, lead iciso strong that 2t 


can be observed in superconducting lead particles. 
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